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— 1 1 
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TECHNICAL MEMORANDUM X-64536 


A CONTROL LAW FOR DOUBLE-GIMBALED CONTROL MOMENT 
GYROS USED FOR SPACE VEHICLE ATTITUDE CONTROL 

SUMMARY 


Space vehicle attitude control, which utilizes control moment gyros 
(CMG's) to develop the necessary control torques, requires the generation of 
CMG gimbal rate commands in such a way that the resulting precessional 
torques on the space vehicle equal the desired control torques; i.e. , no torque 
crosscoupling occurs. Consideration of the combined effect of a pair of double- 
gimbaled CMG's allows the generation of a no-crosscoupling CMG control law 
on the basis of easily understandable kinematic relationships. For the control 
law presented, the only difference between the commanded and the actual con- 
trol torques exerted on the space station is caused by the difference between 
the commanded and the actual gimbal rates. The control law is expanded from 
the application to one CMG pair to the application to three CMG's. Three CMG 
pairs can then be formed and the desired control torque can be split between 
them according to their relative control capability. Skylab-A is used as an 
example for the utilization of the excessive degrees of freedom to better distri- 
bute the CMG angular momentum vectors with respect to each other or their 
gimbal stops without an effect on the total CMG angular momentum; i. e. , with- 
out resulting in a net torque on the space vehicle. The general development 
of the CMG control law assumes arbitrary CMG momentum magnitudes; but it 
is also shown that the expressions can be simplified if it is assumed that all 
CMG angular momentum magnitudes are equal. This simplified version is 
presently used for control of the CMG's on Skylab-A. 


INTRODUCTION 


It is desirable for many space vehicles (especially for an orbiting 
space vehicle like Skylab) [1-4] to have an angular momentum storage device 
on board to accommodate cyclic angular momentum accumulations. This saves 
thruster attitude control fuel and simultaneously allows the reduction of the 
attitude error. Often three double- gimbaled control moment gyros (CMG's) 
are used. Then the need arises to command six gimbal angle rates to 
create a control torque on the vehicle which matches the commanded torque. 


CMG control laws comtemplated in the past such as the cross product steering 
law [1,2] resulted in crosscoupling; i.e. the actual torque deviated from the 
commanded torque in magnitude and direction, even when ideal 1 CMG's were 
assumed. This report shows that the crosscoupling can be eliminated from the 
control of the CMG's by a law which considers the CMG's always in pairs, 
under the assumption that the CMG's are ideal. 

For convenience, the control law is broken down into a steering law 
(which is the control law proper, and the only one to result in a net control 
torque on the vehicle) and two rotation laws. The conventions and the 
nomenclature of Skylab-A will be used throughout the development. The fact 
that failure of a single CMG necessitates two-CMG operations has also been 
kept in mind throughout the development. 

DEVELOPMENT OF A NO-CROSSCOUPLING STEERING 
LAW FOR A CMG PAIR 


Attitude control of a space vehicle is always achieved by application 
of a control torque T. on the vehicle. For a CMG system with a total angular 

momentum 11^, the relationship holds 



( 1 ) 


where is the change rate of with respect to inertial space. The problem 
is therefore how to effect the desired CMG angular momentum change rate H^,. 

The assumption is made that each CMG has a fixed, though arbitrary, 
angular momentum magnitude, generally different from the magnitudes of the 
other CMG's. Elimination of the crosscoupling in the CMG steering law requires 
that the actual angular momentum change is equal to the desired momentum 
change under the assumption that the commanded and the actual gimbal rates 
are equal. While one CMG cannot satisfy this condition, it is relatively easy 
for a pair of CMG's. This will be shown on pair A (CMG's 1 and 2) as an 
example; in the next section the steering law will be expanded to the other 
possible pairings. Hj and are the angular momentum vectors of CMG's l 
and 2, with the initial positions indicated by the subscript I and the final by the 


1. A control moment gyro which has no gimbal inertia, whose angular 

momentum magnitude and direction are known exactly, and which follows 
the commanded gimbal rates exactly. 

2 


subscript F. The desired momentum change is shown in Figure 1 as a momen- 
tum difference AH between the initial pair sum H and the final pair sum 


H . The angular momentum vectors are all shown lying in the same plane 
SA F 

as H_ and II This was done for clarity, but it would generally not be 

dA 1 oA r 

the case. On the other hand, H_, . is not disturbed by a rotation about itself 

— oA 


and a rotation is therefore permissible for the development of the momentum 
change (it is indicative of the fact that one degree of freedom remains, which 
will be treated later) . 



Figure 1. Momentum change. 


The change of H will be broken down into a rotation 0 of H (a 
— SA oA 

rotation of the CMG vector pair as a unit) and into a change in magnitude of 

H by a change of the angles and a 2 (scissoring action of the momentum 
—SA 

vectors with respect to each other). Of course, both motions occur simulta- 
neously. It might be of interest to note that the sum H has not only an upper 

limit (Hj + H 2 ) when the angular momentum vectors are parallel, but also a 
lower limit ( |H t - H 2 |) when the vectors are antiparallel. The latter becomes 
important for two-CMG operation. 


3 






Before the angular velocities for pair rotation and scissoring are 
developed, it is convenient to define the following quantities (a bar below a 
letter indicates a vector; a quantity without a bar indicates either a scalar 
or a vector magnitude) ; 


— SA ' 

s Hj + H 2 


[Nms] 

pair sum 

(2) 

h da 

= Hi . H 2 


[ (Nms) 2 ] 

pair dot product 

(3) 

—PA 

= HjxH2 


[ (Nms) 2 ] 

pair cross product 

(4) 

T 

-t-CA 

II 

> 


[Nm] 

torque command (equiva- 
lent to desired momentum 
change 2 ) 

(5) 

s 

Ssa /h S a 



unit vector along pair sum 

(6) 

^p" 

ap A / H PA 



unit vector perpendicular 
to both Hj and H 2 

(7) 

^x s 

<SSA X ^CA^'SSA 

xTc A i 

unit vector perpendicular 

to both H and T 

— SA — CA 

(8) 






t cap 

% ' — CA 

^SA 

[Nm] 

component of along 

pair sum 

(9) 

t cax 



[Nm] 

component of T 

LA 

(10) 


perpendicular to pair sum 


i = 1,2,3 ) (11) 

j = 1,2,3 
k = A, B, C 

m = i,3 ) (14) 

2. Note that a positive torque command for the CMG's results in a negative 
torque (reaction) on the vehicle. 


( 12 ) 

subscripts used 
throughout (13) 
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The pair rotation will be proportional to T CAX and the necessary 
angular velocity command co xA will be along the unit vector u^. The effective- 


ness of <jo is proportional to H and we get 
XA uA 


-XA T CAXV H SA 


(15) 


or 


— XA SA X — CA^ H SA 


(16) 


Figure 1 will be used as an aid in the development of the angular rate 

command co and co needed for scissoring. They will be proportional 
lr A 1 FAZ 

to the component T = H . The angular velocity for scissoring will be 
CA Jr oA 


co 

— PA1 


O'! U 

-p 


and 


LO 


PA 2 


- cc, u 


“P 


(17) 


(18) 


The following two equations hold 
H 1 ca l + H 2 ca 2 = H gA 

+ H 2 sa 2 = 0 

Differentiation yields 


— HjSQ! | 

-H 2 sq; 2 


“l 


1 

$ 

Hjccuj 

H 2 ca 2 


«2_ 


0 


“1 = H 2 c a2 H gA /H pA 
= -H 1 cc< 1 fl SA /H pA 


(19) 

( 20 ) 


( 21 ) 


( 22 ) 

(23) 
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with 


H pA = HiH^sajca-i - co^sai) 


= - Q!i) 

= |Hjl x Hzl 

With the relationships 


(24) 


H 


SA 


a- ug = 

(H , 2 A Hj • 

(25) 


(H 2 + H, . H2)/H sa 

(26) 

t cap 


(27) 


and equations (22) and (23), equations (17) and (18) become 
— PAl = * H DA ) ( -SsA ' Ica'WVSa > 

- (H * + “da' <2sa ' Ica>!V (h sa h pa ! > 


Cl) 

PA 2 


(28) 

(29) 


Both CMG's participate in the pair rotation through [equation (16)] and in 

the scissoring through u) pAl and I e( l uations ( 28 ) and (29) ] such that the 

angular velocity commands are 

(30) 


St = — xa + —pa i - & 


PAl 


W2 — XA + —PA 2 


(31) 


The angular velocity of the vehicle must be subtracted since — PAl’ anc * 

a; . are with respect to ine 
—PA 2 

with respect to the vehicle) . 


a; are with respect to inertial space (otherwise H_ # ^T) , but and cl^ are 
PA2 
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EXPANSION OF NO-CROSSCOUPLING STEERING 
LAW TO THREE CMG's 


When three CMG's are operative, they can be paired three ways: H 1 
and H 2 form pair A, H 2 and Hj form pair B, and H 3 and H t form paii; C. Each 
CMG participates in two pairings and the resulting angular velocities must be 
added. Basically each pair can produce the commanded torque and a means 
must be found to split the total command into individual pair commands in such 
a way that the pair capabilities are considered. Equation (16) shows that co 

-XA 

is proportional to 1 /H and equation (28) or (29) shows that the w ,'s are 

bA RAi 

proportional to 1/H . . Since H goes to zero when H reaches its minimum 

Jr A irA oA 

(zero for H t = H 2 ) , the splitting or prorating will be done with a function of 
H (prorating must be identical for co and co . of the same pair) . While 

Jr A K Jr K.1 

prorating with H directly will make the cc's insensitive to H when H 

-tK. JrK JrK 

goes to zero, it is desirable to have the angular velocity commands go to zero too, 
so that the case of one-CMG-out (i. e. , failure of a single CMG) can be accepted 
without modification. Prorating is therefore done with H 2 

JrK 

Before the variable gains used for prorating (splitting) of the torque 
command are developed, it is convenient to add the following definitions. 


SSA 

- Mi 

+ M2 

[Nms] 

(32) 


- M2 

+ M3 

[Nms] 

(33) 

—sc 

= H 3 

+ Mi 

[Nms] 

(34) 

Jix 

= Mi 

+ M2 + M3 

[Nms] 

(35) 

V 

- Mi 

• M2 

[ (NMS) 2 ] 

(36) 

h db 

- M2 

• M3 

[(Nms) 2 ] 

(37) 

h dc 

= M2 

• Mi 

[ (Nms) 2 ] 

(38) 

—PA 

= Mi 

x M2 

[ (Nms) 2 ] 

(39) 
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^pb 

= H 2 XH 3 

[ (Nms) 2 ] 

(40) 

—PC 

= H 3 x li 

[ (Nms) 2 ] 

(41) 

ZH p 2 

‘ h pa + h pb 2 + V 

[ (Nms) 4 ] 

(42) 

With the above definitions and the preceding discussion, 
gains become 

the variable prorating 

g a = 

h pa 2/sh p 


(43) 

g b - 

h pb 2/zh p 


(44) 

G c = 

h pc 2/sh p 


(45) 

and the torque commands become 



T 

-CA 

= g a t 

A — C 


(46) 

T 

— €B 

= G T 
B-C 


(47) 

T 

— CC 

= G T 

c — C 


(48) 


The angular velocity commands for pair A from equations (16), (28), and (29) 


are now 


w, . = (1/H 2 ) (H A x T A ) 

— XA v SA ' v — SA — CA 


(G A /H 2 ) (H x T ) 
v A SA ' SA — C' 


- L H i* /< V a, p>J<S B A’'V 

2p A1 - [<«’ + h da' ( Ssa • Ica' /( V h pa 2) ] £b 

= [g a (H, 2 + H ) (H A • T)/ (H 2 H 2 )1 H 
L A v 2 DA' "■ — SA — C ' SA PA 'J — I 

= ["(H , 2 + H )(H 0 • T _) / (H 2 SH 2 )1h 

|_ 2 DA' v — SA — C" v SA p'J— PA 
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The angular velocity commands for all pairs are then (pair A commands are 
repeated for completeness) 


— XA = ^P A Z// ^ H S A ! ^ H p 2 ^] SA X — 

%b = [ h pb 2/(h sb ! sh p>] %b x V 
Sxc - [ h pc , /< h sc ,zh p»J<^o x ^ 

= [(H 2 2 + H 
= [(Ha 2 + H 
= [(H! 2 + H 

^p A2 - + h da>'V • ic ,/(H s A !2H P>]a 

— PB3 = -[ (HZ * H DB ,( SsB • Ic ,/(H SB ZSH P>' 

—pci - -[ (H » 2 + h dc )( Ssc • Ic )/(H sc 2 “p* 


cu 

— PA1 


U> 

— PB2 


CO 

—PC 3 


D A ,( SsA • Ic> /(H SA 2£H P>]iip A 
DB* %SB ’ 2<;^ tH SB 2 2H p 2, j— PB 

dc HHsc • Ic>/< H sc 22 H P 2) ]fipc 


PA 


H 

— PB 


^PC 


(52) 

(53) 

(54) 

(55) 

(56) 

(57) 

(58) 

(59) 

(60) 


The CMG angular velocity commands resulting from the steering law are 

^XA + — PA1 + — xc + ^PCI - iy (61) 




COo — CO + 00 + 00 + c 0 - (I) 

“XB — PB2 “XA “PA 2 % 


COo = CO + CO +U) +00 - d> 

— 3 — XC —PC 3 — XB — PB3 -‘-V 


(62) 


(63) 


Appendix A shows that equations (52) and (60) can be simplified if the 
angular momenta of the CMG's are equal. 

Appendix B shows that the actual torque T is equal to the commanded 
torque T if the commanded and the actual gimbal rates are equal. 



TRANSFORMATION OF A GENERAL RATE INTO GIMBAL RATES 


The angular velocity commands [equations (61) to (63)] are generally 
not perpendicular to the CMG angular momentum vectors and do not depend 
on the CMG mounting configuration. The gimbal rate commands obviously do 
depend on the mounting orientation of the individual CMG. The CMG mounting 
configuration for Skylab-A is used (Figure 2). This configuration is cyclicly 
permutable and the gimbal rates will be developed for CMG l and then 
permuted for the other two. 

The momentum change of CMG 1 resulting from tne commanded 
velocity uq should also result from 6^ (where it is assumed that the actual and 
the commanded gimbal rates are equal) : 

uq x H x = x Hj (64) 


or 


(uq - 6 t ) x Hj = 0 

Geometric relationships give (s = sin, c = cos) 


+6 1(1) S6 3(1) 


+ i Ml)° 6 3 (l ) 


- 6 


Hi - Hi 


3(1) 


+06 1(1) 06 3(1) 


■ C6 l(l) Sd 3(l) 


-s6 


1 ( 1 ) 


With equations (66) and (67), equation (65) results in 
<1) C6 3(1)K<1> 


-(oo -6 c<5 \s6 . +( oo +6 ,\c5 

V 12 1(1) 3(l)y 1(1) V 13 3 < 1 v 1 


(65) 


( 66 ) 


(67) 


(D S6 3(1) 


= 0 ( 68 ) 
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OUTER GIMBAL 
TORQUER 


83(3) 


Z CMG 
MOUNTING 
PLANE 


Sl(3) 


CMG 3 


OUTER 

GIMBAL 

8,(2) ™ 0UER 


\CMG 1 
1 \\\ 

kOV. 


f / / 
CMG 2 V | 


OUTER \ 
GIMBAL Cl 
TORQUER 


. Y CMG MOUNTING PLANE 


X CMG MOUNTING PLANE 


Figure 2. Control moment gyro orientations. 


(“l3 + V>) 06 l,l) Ca 3 W + (“U - 6 1 (l) Sd 3( 1 )) S6 3(D = ° < 69 ' 

- (“n - ‘KD^an^Kij^W 

- (“12 - { i(i) c6 3(i)) c6 MD c6 3(i) = 0 ' l70) 

Equation (70) yields 


6 1(1) = W) + "l2 C6 3(D 


This result inserted into equation (68) or (69) yields (t = tan) 

*3(1, ' - t6 l,l)(“n° 6 3(l, - “l2 S5 3(l)) - “13 ' 

Putting equations (71) and (72) in matrix form and permuting for CMG’s 2 
and 3 results in 




(73) 


(74) 


( 75 ) 
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The direction cosines of the CMG's are usually available and we define e. to 


— 1 


be a unit vector along the angular momentum vector of the ith CMG with the 
following result (for the Skylab-A CMG configuration, Figure 2): 


*1 = 


11 


12 


13 


+C6 1(1) C6 3(1) 


' C6 1(1) S6 3(1) 


-s6 


1(D 


(76) 


-2 


"21 


22 


"23 


-s<5 


1 ( 2 ) 


+c6 c6 

1 ( 2 ) 3 ( 2 ) 


-° a K2) S6 3(2) 


(77) 


—3 = 


31 


32 


33 


-c6_ s6 
1(3) 3(3) 


-s6 


1(3) 


+C6 1(3) CS 3,3) 


(78) 


Equations (73), (74), and (75) can now be expressed in terms of the e... 
With the additional definitions of ^ 


e 1 = 
13 

1/c6 md 

= 1/Vl - e 13 2 

(79) 

e 1 = 
21 

1,/c6 l (2) 

= 1/n/i - e 21 ! 

(80) 

e 1 = 
32 

1/C5 1(3) 

CM 

CM 

CO 

<D 

1 

S 

II 

(81) 

the result is 




6 1(1) 

= e !3 ,(e il 

c o - e to ) 

12 12 ir 

(82) 

6 1(2) 

= e 2l' <e 22 

co - e to ) 

23 23 22' 

(83) 
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( 84 ) 


5 = e ' (e co -e co ) 

1(3) 32 * 33 31 31 22' 


3(1) 


3(2) 


3(3) 


= (e ') 2 e (e co +e co )-co 
v 13 1 11 11 12 12' > 13 


= (e ') 2 e (e co +e co )-co 

v 21 ' 21 1 22 22 23 23 21 

= (e ') 2 e (e co +e co )-co 

1 32 7 32 1 33 33 31 3l' 32 


(85) 

( 86 ) 
(87) 


where it should be remembered that co. = col(co. , co. co. ) can be any 

, , —l ll i2, i3 J 

angular velocity. ’ 


ROTATION LAWS 


Three of the six degrees of freedom of the CMG configuration are 
used for the generation of a control torque on the vehicle. The remaining 
three are the rotations of the pairs about their sums (and also a rotation of 
all three CMG's together about their total sum which is a linear combination 
of the sum rotations) . All these rotations do not result in a momentum 
change; i. e. , no torque is exerted on the vehicle. The rotations have the 
following form: 


co = 

-RA 

ViV'W 

(88) 

co = 

— RB 

V2sb /h sb> 

(89) 

ii 

o 

3 

£ c ( Ssc /H sc ) 

(90) 

CO = 

-RT 

W H T> 

(91) 


and can be used for some benefit. Contrary to the steering law, there is no 
unique way to determine the epsilons. One possible solution is given which 
proved successful in gimbal stop avoidance and in keeping the vectors well 
separated; the implication being that the CMG's have limited freedom of gimbal 
movement. All epsilons can be used for gimbal stop avoidance (R-subscript) 
but only e , e , and e can be used for a proper distribution (momentum 

vector separation; D-subscript) and they are therefore split into two parts: 


e = e + e 
A RA DA 


( 92 ) 
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(93) 


€ B C RB + € DB 


€_ = + e 

C RC DC 


(94) 


The distribution only applies for the case of three CMG's. No distribution 
is necessary for the two-CMG case (the vectors are already located at their 
proper separation and this separation depends upon the sum which cannot be 
altered) . 


Gimbal Stop Avoidance (Rotation Law) 

Avoidance of the gimbal stops is treated first, and it is referred to as 

the rotation law (R-subscript) in spite of the fact that the distribution law also 

uses rotations about the individual momentum vector sums. CMG pair A is 

used again for the development. Four gimbal angles are affected by e • 

KA 

i. e. , a compromise is necessary, and it is therefore desirable to make 
the sum of the individually desirable rotations: 


€ RA e All + e A31 + € A12 + C A32 


(95) 


A desirable rotation is such that the gimbal angle magnitude is reduced; i. e. , 

6 ... 6 <0. Therefore each component of e . was chosen to be of the 

m(i) m(i) RA 

form 


e 


Ami 



F S 
mi Ami 


(96) 


where K„ is a fixed gain, F . is an odd gimbal angle function, and S A . is a 
R mi Ami 

modified sign function. The F-functions for the inner gimbal angles are 

F li - < 9 ’ 

and for the outer gimbal angles we select 

F 3i - [ R ( 6 3(i,)-’ /4 ] n ' (9: 

The F-functions for the outer gimbal angles had to be modified because the 
center between the stops isat+7r/4 for Skylab-A, which is the example used 
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throughout this report. The multiplier R is the ratio of inner to outer gimbal 
freedom giving the gimbal angles equal weight at their stops. The fifth power 
(n = 5) was found by simulation to be appropriate for three-CMG operation 
where for small gimbal angles the distribution (vector separation) should have 
preference. For two-CMG operation, however, the first power was found to 
be more suitable. 


The need for the S^^-functions arises from the fact that the polarity 


of the gimbal rate depends also on the direction of the pair sum with respect to 
the individual CMG. To establish the S A .-functions, a unit vector along 

F is used with the gimbal rate equations (82) through (87): 
oA 

S AU ~ e i3'[ e u tH SA2 /H SA> ’ e 12 < H SA1 /H S A>] 


= (e 'H /H OA )(e. ,e - e e ) 
v 13 2 SA M 11 22 12 2l' 


S A12 * V[VW*W ‘ e 23< H SA2 /H SA> 


] 


= (e *H /H ) (e e - e e ) 
V 21 r SA ; V 22 13 23 12' 


S A31 = e n« H SAl /H SA> + e i2< H SA2 /H SA> 


- < H SA3 /H SA> 


S A32 < e 2 1 ’ ^ ®21 6 22^ H SA2^ H SA^ + ®23 ^ H SA3 //H SA^ 




- (H SA1 /H SA ) 


T SA' 


= ( h </ H oa) (e 21 ') 2 e 0 , (e 00 e <0 + e 1)0 e 40 ) - e 


21 ' 22 12 23 13 


■»] 


For pair B and C we get (through cyclic permutation): 
S, 


B12 


(e ' H /H ) (e e e e ) 

V 21 3 X SB' V 22 33 - 23 32' 


(99) 


( 100 ) 


= ‘V^[< e l3’>V , ll e 21 + e i2 e 22 ( - e 2 3 ] (101 > 


( 102 ) 


(103) 
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(104) 


S C13 (e 32' H l //H SC ) (e 33 e il e 3i e i3^ 


S B13 (e 32' H 2' /H SB^ (e 33 e 21 ®3i e 23^ 


S CU = (e i3 ,H 3^ H SC^ e ii e 32 " e i2®31^ 

S B32 = **V H SB^ [^ e 2l'^ e 21 (e 22 e 32 + e 23*W " 6 31 
S C33 = ( H l/ H SC ) [ (e 32^ e 32 (e 33 e i3 + 6 3i e il^ “ ®12] 
S B33 = ^ H 2 ,/H SB^ f e 32^ e 32^ e 33 e 23 + e 31 6 21^ ” ^ 2 ] 


'] 


S 


C31 <V H SC ) (e i3^ e !3 <6 1 I e 31 + ®12 e 32^ ®33 


(105) 

(106) 

(107) 

(108) 

(109) 

( 110 ) 


All gimbal angles are affected by a rotation about the total angular momentum 
and we select 


€ T 6 T11 + e T12 + € T13 + £ T31 + € T32 + € T33 


( 111 ) 


Again the selected form is 

= - K_ F . S„ 


'Tmi 


R mi Tmi 


( 112 ) 


The angle functions F . are unchanged. The S ^ terms will be developed 
along the same line as the other S-functions (CMG 1 serves as an example) 

S T11 = e i3'[WV - e i2 (H Tl /H T ) ] 

(e i3' H 2^ H T^ (e ii e 22 “ e i2 e 21 ) 

+ (e i3' H 3^ H T (e ii e 32 " e i2®31^ 


(H SA S A11 + h sc s cii> /h t 


(113) 
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S T31 = ( e i3^ e i3[ e il (H Tl /H T ) + V’WV] " 

(I V*V [ (e i3^ e i3 (e ii e 21 + e i2 e 22^ " e 2s] 

+ (H 3 /H t ) j*( e 13 ') 2 e 13 ( e 11 e 31 + e 12 e 32 ) ~ e 33 J 

(114) 

(115) 

(116) 

(117) 

(118) 


= (H SA S A31 + H SC S C31 )/H T 


Cyclic permutation yields for CMG's 2 and 3 


S T12 (H SB S B12 + H SA S A12^ H T 


S T13 (H SC S C13 + H SB S B13^ //H T 


S T32 (H SB S B32 + H SA S A32^ ,/H T 


S T33 (H SC S C33 + H SB S B33^ //H T 


The S-functions will have an upper limit of S . This provides a 

J-J 

linear range (besides the sign information) which is needed to avoid limit 
cycling otherwise introduced by the sign changes. The value of S allows 

selection of the linear range and the loss in gain (S < 1) can be made up by 
an increase in K . 

rv 

CMG Vector Separation (Distribution Law) 


The bulk of the momentum change occurs quite frequently along a well 
known axis [5], If this is the case, one type of distribution law (D-law) can 
be applied which will separate the CMG vectors by trying to make them contri- 
bute equally (in proportion to their magnitudes) to the angular momentum along 
this axis. 

Let e^ be a unit vector along the bulk of the momentum change (the 
orbit normal for Skylab) ; then the angle between e N and the plane formed by a 
vector pair will be maximized by (for pair A as an example) 
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£ da - k d|[<W - <VV ■ <“ S a /h sa><“sa /h sa>] 

- [<V H N> - <W • <S SA /H SA> ( !isA /H SA , ]|-S N < 119 ) 

The significance of the various terms is illustrated in Figure 3. is a con- 
stant gain and H T is a nominal momentum used for normalization. • Evaluation 
N 

leads to the intermediate step 

£ da - ‘Wfr-3*- (H ' J - h ^i + V /h sa ! ]-% • U20) 


The third term in the brackets shows that it vanishes when the magnitudes of 
the angular momentum vectors are equal. Further evaluation yields 


'DA = <V h nV 


>[« 


H 


h da^i 


< H / + •WSs] • % < 121 > 


with H = H • H [cf. equations (28) and (29)]. It should be noted that 

DA "*"""1 Zt 


SgA * Sp A - <H l * Vil - < H i‘ + 

With this equality and the use of cyclic permutations we get 

e A = (K n / H J D1 2 ) (ILa x H„J • e XT 
DA D N SA SA — PA — N 


DB 


DC 


< VVW* <Ssb x M P b> • % 

< V h n h sc 2) <»sc x Spc> • % 


( 122 ) 


(123) 

(124) 

(125) 


This form results in an isogonal 3 distribution about the total angular 
momentum vector if is along II , if all CMG momentum magnitudes are 

equal, and if there is no rotation law effective (K^ = 0). Otherwise, a com- 
promise results between the tendencies to spread the vectors and to reduce the 
gimbal angles. 


3. A distribution of three control moment gyros of equal momentum magnitude 
which contribute equally to their total momentum. This distribution results 
in equal angles between the individual momentum vectors and the total vector 
and equal angles between the vectors themselves. 
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ACTUAL PROJECTION (Hj /H N )-^ N 

DESIRED PROJECTION [(Hj /H N )*(H SA /H s/ j] (H SA /H SA ) • § N 

NOTE: H SA /H N a e N LIE IN THE PAPER PLaNE 

H, /H n AND ti 2 /H N 00 N0T LIE IN THE PAPER PLANE > 

BUT IN A PLANE WITH H /H M 

“SA N 


Figure 3. Distribution. 
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A positive K R yields a right-handed and a negative yields a left-handed 
configuration (looking down on e , right-handed means a sequence of 
11 ) . Both configurations are stable. The location of the gimbal stops with 
respect to determines which of the two configurations is preferable. 

TOTAL CMG ANGULAR VELOCITY COMMANDS 


The CMG angular velocity commands from the various sources can 
be vectorially added to form the total CMG angular velocity commands: 


CO = C0 + CO +C0 -f CO + O) + CO HK CO - <f) 

— 1 — XA — PA1 — XC —PCI — RA — RC — RT ^ 


■V 


( 126 ) 


CO — to Hh CO + c 0 + CO + CO + CO + (0 - d) 

—2 — XB — PB2 — XA —PA 2 — RB — RA — RT 


( 127 ) 


» 

CO ~ CO + CO + CO + (0 -f CO + 00 + ^ — <f)__ 

—3 — XC —PC 3 — XB — PB3 — RC — RB — RT X V 


( 128 ) 


CONCLUSIONS 


A no-crosscoupling control law for double-gimbaled CMG's can be 
developed using easily understandable kinematic relationships. The control 
law is based on a CMG pair as the smallest unit, which can give a no- 
crosscoupling control law; but the law lends itself to easy expansion to the 
control of any number of CMG's as shown by the expansion to three CMG's. 

The development is based on the restriction that the CMG's have fixed momen- 
tum magnitudes, though the individual magnitudes are not necessarily equal to 
each other. Unequal magnitudes resulted for the two-CMG case in an upper 
as well as a lower limit for the total angular momentum. The lower limit 
approaches zero when the momentum magnitudes become equal, which 
therefore is a desirable characteristic. 
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APPENDIX A 


CONTROL LAW FOR NOMINAL ANGULAR MOMENTUM MAGNITUDE 4 


If the assumption is made that the angular momentum magnitude of 
each of the CMG's is equal to the nominal value H^, simplification and norma- 
lizations can be applied. We have 


“l = h n^i 


—2 = V2 


—3 = H N— 3 


(Al) 

(A2) 

(A3) 


where the e*s are unit vectors along the CMG's angular momentum [equations 
(76) to (78)] whose components usually are available from gimbal resolver 
chains. Using pair A as an example we also have 


— SA ■ + •*> 


= tL T e„ A 
N-SA 


—PA 


H N (— 1 X — 2 ^ 


= h n^pa 


(A4) 


(A5) 


The relation H = H 

1 Lt 

(i = 1,2) 


(H? + 




= H. T leads toa„ 
N 2 

< H N + H N c2a > 
2(H N ca) ! 

2 < h sa /2)2 ■ 


-a 


a and this results in 


(A6) 


4. This simplified version is presently used for control of the CMG's on 
Skylab-A. 
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Consequently we have 


< H i 2 + h da )/h sa 2 = ^ 


(A7) 


Normalized Steering Law 

Introduction of the above normalizations and simplifications yields 

'PA / ( e SA ^ I)] (%A * A) <A8) 


-XA 


with 


A = V h n 


(A9) 


and 


Se J = e 


PA 


+ e FB + 6 


PC 


(A10) 


We now have (co pAl = - o) pA2 = cu pA ) 


Cl) 


PA 


( -SA ' ^ 


ic)/(2^)]£ PA 


(All) 


With 


H = H = H it is now possible that the sum H + H goes to zero 
1 2 N 12 

which will result in co being indeterminate since the cross product goes to 

XA 

zero simultaneously. This can be avoided by the relationship 
(e pA / e SA ) 2 = (s2a) 2 / (2ca) 2 


= s^a 


= 1 - era 


- ‘ - (e SA /2 > 2 


(A 12) 


With equation (A12), co of equation (A8) becomes 

XA 


— XA = [(>-e SA V4)/Se^](e| 


SA x 2c> 


(A13) 
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The pair rate commands for nominal are (pair A commands are repeated 


for completeness) 


= [(l - e SA 74) /Sep] (e SA x 

(A14) 

“XB = [(‘ - e SB /4 ) /26 p] X 2c> 

(A15) 

-xc = [(* " e sc' /4 )' /Ze p] < -sc * 77 

(A 16) 

"p A - [(Ssa • ic)/( 22e p)]£p A 

(A17) 

-PB = [‘-SB ' 77'" ( 226 p)] -PB 

(A18) 

^PC = [<%C ■ ( 2Ze p)] -PC • 

(A 19) 

For the nominal H case (H = H = H = H ) the CMG angular velocity 

1 & o JN 

commands are * 

co = co + co A + UJ - co - d> 

—1 — XA —PA — XC —PC % 

(A20) 

CO — CO + CO + QJ — (jj — cf) 

-2 -XB -PB -XA -PA % 

(A21) 

-3 -XC -PC -XB -PB % 

(A22) 

Normalized Rotation Law 


When normalization is introduced, equations (88) to (91) change into 

— RA = e A ( ^SA /e SA ) 

(A23) 

— RB ~ € B ( -SB^ e SB^ 

(A24) 

— RC " e C ( -SC //e SC^ 

(A25) 

— RT “ e x -T^T^ 

(A26) 
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where 


-T = -1 + -2 + -3 ‘ 


The S-functions are also affected and change into 
S, 


All 


A12 


B12 


B13 


C13 


Cll 


Til 


T12 


T13 


A31 


A 32 


B32 


B33 


C33 


C31 


(e '/e ) (e e - e e ) 

1 13 SA' 1 11 22 12 2i 


(e '/e ) (e „e - e e ) 

' 21 7 SA' v 22 13 23 12' 


(e '/e ) (e e -e e ) 

1 21 7 SB 7 V 22 33 23 32' 


(e '/e ) (e e - e e ) 

' 32 SB' 1 33 21 31 23' 


( e */e ) (e 3 -e e ) 

V 32 7 SC' V 33 11 31 13 


Ce '/e ) (e e -e e ) 

' 13 7 SC' 1 11 32 12 3l' 


(e SA S AU + e SC S C 11)/e T 


(e S + e S )/e 
V SB B12 SA A12 ;/ T 


(e S + e S )/e 
K SC C13 SC C13 ;/ T 


r(e ! ) e (e e + e e ) 

'A O ' A O ' i H Ol 4 Q OO' 


13 13 11 21 12 22 
L'~21 ' 3 21^ e 22 e i2 + e 23 e i3^ 

[* e 2l'^ e 21 (e 22 e 32 + e 23 e 33^ 

Re ') e (e e + e e ) 

32 ' 32' 33 23 31 2l' 

Re ') e (e e + e e ) 

L 32 ' 32 ' 33 13 31 ll' 

Re ') e (e e + e e ) 

L' 13 ' 13 ' 11 31 12 32' 


e 23] /e SA 


e il] /e SA 
e 3l] //e SB 
e 22] /e SB 
e i2]^ e SC 
e 33] /e SC 


(A27) 
(A28) 
(A29) 
(A30) 
(A31) 
(A32) 
(A33) 
(A 34) 
(A35) 
(A36) 
(A37) 
(A38) 
(A39) 
(A40) 
(A41) 
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(A42) 


S T31 (e SA S A31 + e SC S C31^ e T 


S = (e S + e S )/e 
T32 v SB B32 SA A32 ;/ T 


(A43) 


S = (e S + e S ) /e 

T33 1 SC C33 SB B33 7 T 


(A44) 


Ail other equations stay the same [equations (95) to (98) , (lll)and (112) ] 


Normalized Distribution Law 

Equation (120) of the distribution law allows easy change for the case 

that H. = H which results in 
x N 


DA 

- k d% 

-v 

• % 

(A45) 

DB 

■ K i>% 

- V 

' % 

(A46) 

DC 

II 

|o> 

CO 

- V 

• % • 

(A47) 


27 


APPENDIX B 


PROOF OF TORQUE EQUIVALENCE 


In the following development the proof is given that the actual torque T? 
is equal to the commanded torque under the assumption that the commanded 
and the actual gimbal rates are equal. 

Several identities are needed for the development and are given first. 


(Vj x V 2 ) x V 3 = (Vj • V 3 )V 2 - (V 2 
where the V's are arbitrary vectors. 

H PA ■ H / H 2 s2 <“2 - V 

= H ! !H 2 - H 1 ! H 2 C2< “2 - V 


= “l H 2 - H DA 


V,)V, 


(Bl) 


(B2) 


The torque applied to the CMG's is (the opposite of this torque is the 
torque on the vehicle) 

T = + _£ y ) x H t + (o> 2 + ± y ) x H 2 + (co 3 + x H 3 . (B3) 


With equations (61) and (63) we get 


-XA + 

CO 4- 

10 + 

— PCl^ x 

H 

— PA1 

— xc 

— 1 

+ XB 

+ CO 

+ CO 

+ — PA2^ 

X H n 

— PB2 

XA 

—2 

+ texc 

4* CO 

+ cu 

+ — PB3^ 

x H 

—PC 3 

— XB 

—3 


(B4) 


or 


28 



I = Hxa x ®1 + S2» + — xb x <£ 2 + h 3 ) + "xc x <Ss + £i> 


+ Sip A1 X Sj + “p A2 x —2 + Hp B2 X S 2 + — PB3 X “3 
+ — PC3 X —3 + “PC 1 x — 1 ■ 

Further development with equations (52) to (60) and (Bl) yields 


1 = L H pa 2/ ( h sa 


2 SH r 2 'M (H^ a x T> x H 


— SA — C 


■SA 


) 


+ [ h pb / ( h sb sh p)] <Ssb x V x Ssb 
+ [ h pc 2/ ( h sc sh p)] ‘Ss C x Ic )x £sc 


'( H 2 + H Da) %A ' Ic (/ ( H SA ^p)] <£l X £ 2 > X £l 


[(”< 


+ h da) <«sa • V/ ( H SA ZH " 


): 


® 2 xa 1 )x a 2 


(B5) 


( H 3 + H Db) <£SB ' So»/ ( H SB 2H ;)] X S*> X £ 2 

+ [ («2 + H Db) «SSB ' 2c’/ ( H SB ZH p)] <«3 X £ 2 » X £3 

+ [( H 1 2 + h dc) ( £sc ' 3c' 7 ( h sc sh p)] <£3 x V x £3 

+ [(»3 + h dc) <£sc • 3c J/ ( H sc sh p)] <£i x £s» x £1 (B6) 

3 = [ h pa 2/ ( h sa sh p)] [ h sa 2 3c - ( £sa ' 3^] 

+ [ h pb 2/ ( h sb 2sh p)] [ h sb 2 3c - <£sb • 3c)Hsb] 

2^u 2V1 r H _ ffl • T ) 1 1 1 

-c — sc] 


+ [ h pc 2/ ( h sc sh p)][ h sc 2 3 c - ( £sc 


( , 2 +B nd)®«-2c J 


m 3 + h db) ( £sb • 3c) 


+ H. 


s) 


(H 


DBI ‘-SB -C' 


/ H 2 2H 2 \ 

'] 

\ SA P /_ 


(k 2 SH 

] 

\ SA P^J 


( h sb sh p)] [' 

(H 2 SH 2 V 

[] 

\ SB Pyj 


/h_ 2 sh 

["] 

V sc p/_ 

L 

(h 2 sh 2N \1 [1 
V sc P/J _ 


h i 2 £ 2 - H DA—1 


2—1 DA-2 


-3 - H DB 




[ 3 2 £2 


[C 

[(") 

_( H 1 2 + h dc) { £sc • 3 c) 

K* h dc) <£sc • 3c)/ ( h sc eh p)] [ h i-3 - h dc£i] (b7) 


DB— 3 
3-1 ^DC-3 
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T = [Yh a 2 +H 2 + H 2 VsH„1t_ 
\\ PA PB PC/ Pj-C 


2 SH 2 


T* \ u 

SA — C — SA 


- [V^i 

- [ h pb V ( h 

- [ h pc V ( H sc sh p)] l “SC ' VS 


„ 2 SH 2 


1 
)] 


<Ssb ' V^b 


-sc 


(»2 “l^— 2 - H 2 H DaS, + H DA H l ! a 2 - VSl 


+ H/H 2 ^ - H/H da H 2 ♦ H da H 2 ^ - H da >H 2 ) (J!sa • T, 
+ ( H 3 H 2 ! ii 3 - H 3 H DB^2 + “dB^S ' "mk 

+ H 2 H 3ii 2 - H 2 H DBfi 3 + H DB H 3 a 2 ' “dB^) '“sb' $ 

+ ( H 1 H 3% - H ,\c^ + H DC H 3fi! - H DC^3 
+ H S H 1% ~ H 3 H D(A + H DC H 1^3 

-V^) ( 2sc ' ( H sc ZH p ! ) 

Equations (42) and (B2) allow the reduction of equation (B8) to 


- [ H PA !/ ( H SA ZH p)] <JisA • — C^— SA- 

- [ h pb 2/ ( h sb Zz:h p)] 'Ssb'Ic'Ssb 

- [ H PC !/ ( H sc ^p)] ( Ssc • Ic'Ssc 

- [V«»- V / ( H SB S “p)]a S A 

+ [VSsB-Ic^^SB^pijiisB 
+ Kc'Ssc-Ic^iv^p)]^ 


^("sA^p) 

^(“sB^p) 

(B8) 

(B9) 
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showing that 



which was to be proven. 


(BIO) 
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